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Abstract 


We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec Integral Field Unit 
(IFU) spectroscopy of the young Galactic supernova remnant Cassiopeia A (Cas A) to probe the physical conditions 
for molecular CO formation and destruction in supernova ejecta. We obtained the data as part of a JWST survey of 
Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO 
on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer 
layers than the Ar ejecta, which indicates the re-formation of CO molecules behind the reverse shock. NIRSpec-IFU 
spectra (3—5.5 um) were obtained toward two representative knots in the NE and S fields that show very different 
nucleosynthesis characteristics. Both regions are dominated by the bright fundamental rovibrational band of CO in 
the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si Ix], [Ca IV], 
[Ca V], and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity 
space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal 
patterns of rovibrational lines with pseudocontinuum underneath, which is attributed to the high-velocity widths of 
CO lines. Our results with LTE modeling of CO emission indicate a temperature of ~1080 K and provide unique 
insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae and have strong 
ramifications for modeling dust formation that is led by CO cooling in the early Universe. 


Unified Astronomy Thesaurus concepts: Supernova remnants (1667); Metal-containing molecules (2258); CO line 
emission (262); Supernovae (1668); James Webb Space Telescope (2291); Ejecta (453); Astrochemistry (75); Dust 
formation (2269) 


1. Introduction 


Original content from this work may be used under the terms The 1 u^ of dust : high- laxi 

BY of the Creative Commons Attribution 4.0 licence. Any further . oon AMOUS: O ed PELIN Some lgh-z galaxies 
distribution of this work must maintain attribution to the author(s) and the title imply that dust formed in the early Universe. However, 
of the work, journal citation and DOI. intermediate-mass stars, thought to produce most interstellar 
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dust when on the asymptotic giant branch in present-day 
galaxies (Laporte et al. 2017), would not have evolved to the 
dust-producing stage in high-z galaxies. In contrast, core- 
collapse supernovae (ccSNe) from high-mass stars occur just 
several million years after their progenitors are born and 
have also been suggested to be molecular factories in the 
early Universe (Cherchneff & Lilly 2008). Molecules, e.g., 
carbon monoxide (CO), are the signature of the onset of dust 
formation since they are one of the ejecta's most powerful 
coolants (Cherchneff & Lilly 2008) at temperatures where dust 
can form. 

However, whether SNe are a significant (if not dominant) 
source of dust in the early Universe has been and continues to 
be debated (e.g., Nozawa et al. 2006; Cherchneff & Dwek 
2009). The dust masses observed in ccSNe within a few years 
after their explosions are less than 0.01 Mo, a value that is far 
too small to explain the amount of dust observed in the early 
Universe (Kotak et al. 2009; Gall et al. 2011; Tinyanont et al. 
2019). In contrast, Herschel and Spitzer observations of a few 
young supernova remnants (SNRs), including SN 1987A 
(20 yr; Matsuura et al. 2011) and Cas A (~350 yr), have dust 
masses of 0.1-1 M. (Rho et al. 2018b; Millard et al. 2021; 
Matsuura et al. 2015; Chawner et al. 2019; De Looze et al. 
2017, and references therein), which is in agreement with dust 
formation models (Nozawa et al. 2003; Todini & Ferrara 2001; 
Sluder et al. 2018) that suggest that SNe could be major dust 
factories at high-z galaxies (Dwek & Cherchneff 2011). 

The cause of discrepancies in the measured dust masses from 
early to later phases and the timescale of dust formation are 
under debate. Recently, Niculescu-Duvaz et al. (2022) and 
Shahbandeh et al. (2023) suggested that a dust mass grows with 
time and most dust forms at times >3 yr after the explosion, 
while Dwek et al. (2019) suggest that the dust forms early in 
optically thick clumps and only a fraction of its IR emission is 
detected due to high IR opacity. 

Dust may undergo complete or partial destruction following 
its initial formation, which depends on shock velocity and grain 
size (Slavin et al. 2020), potentially altering the expected dust 
grain size distribution in the process. When the forward shock of 
a supernova (SN) accumulates sufficient pressure at the SN shell, 
a second shock (called a reverse shock) develops in the interior 
(Chevalier 1977; Borkowski & Shull 1990). The dust (or ejecta) 
destruction by the passage through the reverse shock of ccSNe 
depends on grain composition, grain size distribution, and the 
shock properties (Priestley et al. 2021; Nozawa et al. 2007; 
Kirchschlager et al. 2019, 2023). Because the dust formed in 
ccSNe includes sufficiently large grains (0.1-0.5 um), a 
significant fraction of the grains can survive (1096—2096 for 
silicate dust and 3096—5046 for carbon dust; Slavin et al. 2020). 

Cas A is one of the youngest (~350 yr) and closest ejecta- 
dominated SNRs and has been observed across all wave- 
lengths, including optical (Fesen et al. 2006), X-ray (Hwang 
et al. 2004), radio (DeLaney et al. 2014), and infrared (Isensee 
et al. 2010, 2012; Smith et al. 2007). The progenitor mass of 
Cas A is still uncertain; it has been suggested it was a Wolf- 
Rayet star with a mass of 15-25 Mo (Fesen 2001; Young et al. 
2006). Spectra of a light echo from the SN that produced Cas A 
showed a Type IIb SN from the collapse of the helium core of a 
red supergiant that had lost most of its hydrogen envelope 
(Krause et al. 2008). Koo et al. (2020) alternatively suggest a 
blue supergiant precursor with a thin hydrogen envelope or a 
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yellow supergiant (the progenitor mass is probably «15 Mọ, 
depending on pre-SN mass loss). 

Cas A is the best case study in the local Universe to 
understand dust formation in SN ejecta and shock-processing 
of freshly formed SN dust (De Looze et al. 2017). Spitzer 
infrared spectral mapping of Cas A confirmed that molecules 
and dust are present in the ejecta (Rho et al. 2008). Cas A also 
shows strong polarization fractions (~20%) in the far-infrared, 
implying that the dust grains are large (Rho et al. 2023). CO 
was detected for the first time from the Palomar and Spitzer 
near-infrared (NIR) images (Rho et al. 2009, 2012). CO 
fundamental band features were detected in low-resolution 
AKARI spectra, showing that astrochemical processes and 
molecule formation continue to the stage of young SNRs (Rho 
et al. 2009, 2012). 

Although CO molecules themselves are not on the 
immediate chemical path to dust formation (Sarangi & 
Cherchneff 2015), it is known that the detection of CO is an 
indication of molecular cooling and chemistry in the ejecta, 
leading to condensation of ejecta dust in later (from a few tens 
to a few hundred days after the explosion) epochs (Sarangi 
et al. 2018; Rho et al. 2018b, 2021). Although CO can form in 
the ejecta in early phases (~100 days after SN explosion; 
Sarangi & Cherchneff 2013), the detection of later epoch CO 
has been associated with re-formation in the post-shock region 
(Biscaro & Cherchneff 2014), highlighting the complex and 
competing processes of molecular formation and destruction. 

CO molecular lines have been observed in the late-time 
NIR spectroscopy of a number of SNe, e.g., two Type IIP 
SN1987A (Spyromilio et al. 1988) and SN 2017eaw (Rho et al. 
2018a; Tinyanont et al. 2019) and a Type Ic SN202001 (Rho 
et al. 2021). By comparing the production of molecules, ejecta, 
and dust at these early times to the molecular gas in Cas A, we 
can constrain both production and destruction of dust and 
molecules in SN explosions. 

In this Letter, we present new observations of Cas A using 
JWST, focusing on NIRCam and MIRI imaging of the entire 
SNR and NIRSpec Integral Field Unit (IFU) observations of 
two selected filaments. JWST reveals that the CO-emitting 
regions mostly coincide spatially with the ejecta-dominated 
areas; however, the ratio between the CO and Ar ejecta varies 
across the SNR, and fine-scale differences in position and 
morphology exist between the CO emission and the ejecta. The 
NIRSpec spectra are dominated by complex CO bands and 
show a unique shape of the CO fundamental band. 


2. Observations 


We present JWST infrared observations of Cas A using 
NIRCam (Rieke et al. 2023), MIRI (Wright et al. 2023), and 
NIRSpec (Jakobsen et al. 2022). These observations were 
obtained as part of a Cycle 1 survey program on CasA 
(Prog. ID of 1947), described in Milisavljevic et al. (2024). The 
NIRCam images were taken on 2022 December 5-6, using the 
F356W (3.140—3.980 um) and F444W  (3.880—4.986 um) 
filters. The field of view (FOV) of NIRCam is 2/2 x 4.4 with 
the pixel scale of 07065. In contrast, Spitzer images at these 
wavelengths have a spatial resolution of 3" (Fazio et al. 2004). 
We required a 2 x 3 mosaic to cover the full extent of Cas A. 
The NIRCam data required 1.9 hr on source exposure time for 
F345W and F444W images. 

MIRI imaging using the F770W filter was taken multiple 
times between 2022 August 4 and October 26. The MIRI 
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Figure 1. JWST F356W (3.56 jum; top left), FA44W (4.44 um; top right) images, CO fundamental image (produced from F444W image after subtracting the 
synchrotron emission using F356W image; bottom left), and equivalent CO image using Spitzer data (bottom right). The JWST images show far more complexity. The 
CO emission is from the CO fundamental bands and overlaps spatially with the knotty ejecta structures. Faint, extended, diffuse F444W emission seen in projection 
toward the center, hinted at with the Spitzer data, is more noticeable. The image is centered on R.A. 23°23™26565 and decl. +58°49! 14/98 (J2000) with an FOV of 


6/4 x 6/4. The units in the color bar are in MJy sr- !. 


imaging FOV is 1/23 x 1788 with the pixel scale of 0711, 
requiring a 5 x 3 map to cover the SNR. The JWST mosaics 
were astrometrically aligned using the JWST Alignment Tool 
(Rest et al. 2023). The resulting aligned mosaics are shown in 
color in Figures 1 and 2. 

Each NIRSpec-IFU position covers a region 3"7 x 4" in size 
with 0/1 pixels. Four IFU positions were observed for the 


program, as described in Milisavljevic et al. (2024). Here we focus 
on two positions that were observed, one on 2022 November 5 
and the other on 2022 December 10. The position in the north 
(R.A. = 3507873, decl. = 587842) targeted an ejecta knot, and the 
position in the south (R.A. — 3507875, decl. — 587790) targeted a 
particularly bright region of CO emission. These positions (P1 and 
P3 in Milisavljevic et al. 2024) are shown in Figure 3. The other 
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Figure 2. JWST three-color mosaicked images of synchrotron emission (F356W in blue), CO (synchrotron-subtracted FA44W in green), and Ar ejecta (F770W in 
red). North is up, and east is to the left. The detailed structures of the three images are noticeably different from each other. The synchrotron emission (blue) shows 


smooth structures and is dominant outside the main shell. 


two regions are reported by I. De Looze (2024, in preparation) and 
Milisavljevic et al. ( ). We extracted the JWST NIRSpec-IFU 
spectra from each pixel, and the total ~1368 spectra are from 
36 x 38 detector pixels. However, note that the IFU had 53 x 55 
image size after World Coordinate System rotation. The effective 
integration and total exposure times are 145 and 583 s for the N 
field and 218 and 875 s for the S field, respectively. More details 
of the JWST observations of Cas A are described in Milisavljevic 
et al. ( ). 


3. Results 
3.1. Imaging 


In Figure 1, we present the mosaicked NIRCam images. The 
top two panels show the separate JWST NIRCam mosaics from 
F356W and F444W. Much of the F356W emission appears 
relatively diffuse, similar to synchrotron radiation seen in the 
radio (Anderson et al. ), K-continuum (Rho et al. ) 
and Spitzer 3.6 um (Rho et al. ) images. However, we note 
some knotty structures, likely ejecta, in the F356W image of 
the northern shell, which are visibly more compact and brighter 
than the synchrotron-emitting regions. The NIRSpec spectra 
also show a continuum from synchrotron emission, but at some 
locations the F356W band includes a [Ca IV] 3.2 um line (see 
Figure 4 and Section for details). Thus, while there may 


be isolated, localized contributions from ejecta emission, the 
F356W image mainly shows synchrotron emission. 

The F444W image shows three main structures: (1) knotty 
structures from ejecta, (2) the forward-shocked material outside 
the bright shell, which shows a much smoother distribution, 
and (3) weak diffuse structures in the interior (see Figure 1). 
The NIRSpec-IFU spectra show CO emission is dominant 
(29946) in compact knots and filaments, and the F444W filter 
also includes ejecta lines of [SiIX] at 3.94 jum, [Ca V] at 
4.16 um, [Mg IV] at 4.49 um, [Ar VI] at 4.53 um, and [K m] at 
4.6 um (see Section for details). The F444W filter also 
includes synchrotron emission since the forward-shocked 
material is known to be from synchrotron emission, as seen 
in radio images (Anderson et al. ). 

To remove the synchrotron component to first order and 
obtain a CO image, we scaled the F356W image to the F444W 
image assuming a synchrotron spectral index a = —0.68 (Rho 
et al. ) where log S, xa log v and subtracted it from the 
F444W image, which resulted in the lower left panel in 
Figure 1. This image is dominated by the inner ring of emission 
and should be due primarily to fundamental CO band emission, 
as previously suggested using Spitzer IRAC band 4 (at 4.5 jum) 
imagery and AKARI spectra (Rho et al. ); this will be 
confirmed and discussed further using NIRSpec in Section 
The fourth panel in Figure 1 shows the Spitzer CO image for 
comparison (Rho et al. 2012). Individual CO structures in the 
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Figure 3. Zoomed-in images of Figure 2 on the northern (top) and southern (bottom) shells. The NIRSpec-IFU FOVs are marked as black squares on the JWST three- 
color images. The arrows point to the filaments showing CO excess emission (in green, marked as A, B, and D) and Ar ejecta excess (in red, marked as C). The slit 
(box E, in cyan) is cut through the southern shell where a radial profile is obtained in Figure 7. 


1x10* 


6x10*E T zt T TT - T E 

E E gOET zz j 

E 3 = OG Bü j 

E o [s OO zm s 
5x10*F- Es i Ar! TE 4 

E ] ] Tel i j 

~ E Il 1 Fd I j 
a E I | Fg | J 
u 4x10 E f wo Gi i E 
3 F | | ug | J 
— E | | FSI J 
2 E f | l: iH f j 
© 3x10*E- I l Ha | E 
a3 E | | te iy | J 
= E l \ eae l J 
& E ] ] pde I j 
, exiotE l ] D l 4 
8 E | | i | j 
& E | i j l j 
3 E Hoan 3 
n C | | I: -j 
E ] ] H 1 

E | | ] J 

E | : J 

E MM. bi 4 


3.0 3.5 4.0 4.5 5.0 
Wavelength (um) 

Figure 4. JWST NIRSpec-IFU spectra of Cas A toward two representative 
knots (0/1 x 0.1 region for each) in the northern IFU field (R.A. = 
350?8731, decl. = 58°8422 in blue) and at the southern IFU field (R.A. = 
350?8764, decl. = 5877905 in green). The spectra show a dominant CO 
fundamental band with a few high-ionization ejecta lines indicated. The two 
CO bumps arise from the R (the shorter wavelength bump) and P (the longer 
wavelength bump) branches of CO. The CO emission is stronger in the 
southern position, as expected from the imagery presented above. The ejecta 
lines are from [CaIV] at 3.2 um, [Si IX] at 3.94 jum, [Ca V] at 4.16 jum, 
[Mg IV] at 4.49 jum, [Ar VI] at 4.53 um, and [K Il] at 4.6 jum. 


Spitzer image are resolved into numerous knots in the JWST 
images. 

With this information in hand, we show the resulting aligned 
mosaics in color in Figure 2. Synchrotron (F356W) is shown in 
blue, fundamental CO emission (F444W after synchrotron 
subtraction) is shown in green, and the Ar ejecta emission 
(from [Ar II] at 6.98 uum) is shown in red. The Ar ejecta image 
is from the MIRI F770W data,k which cover 6.58—-8.687 um. 
The Spitzer spectrum shows that the Ar lines are dominant 
within these wavelengths (Rho et al. 2008; Smith et al. 2009). 
The synchrotron emission arises mainly from Fermi accelera- 
tion in the forward-shock region as it sweeps up and encounters 
interstellar medium/circumstellar medium (Rho et al. 2003). 
The three-color mosaicked images contrast synchrotron emis- 
sion, CO molecules, and Ar ejecta emission in great detail on 
all scales all the way down to individual ejecta knots. While the 
ejecta and CO emissions both arise primarily in the bright ring 
of emission caused by the reverse shock, their spatial 
distributions are quite different from each other. Here we 
highlight two regions: the so-called *Ne-moon" regions (Figure 
8 of Smith et al. 2009; Ennis et al. 2006) where the Ne 
emission is brighter than Ar emission; these are the two regions 
that appear in green and white in Figure 2. 

We find that the spatial distribution of Ar (eastern region in 
red of Figure 2) and CO regions (northern and southern regions 
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Figure 5. (Left) Two representative NIRSpec-IFU spectra from the northern shell (in orange) and southern shell (green) show dominant features from CO fundamental 
bands with a few high-ionization ejecta lines, in comparison with Orion (blue; Peeters et al. 2024). The northern spectra show two velocity components of —2070 and 
—1270 km s^, while the southern spectra show one component of +1060 km s^ !. The ~3400 km s^ ! difference between the blue- and redshifted motion of Cas A 
has been removed using the rest wavelengths of the ejecta lines. (Right) An enlarged section from the left panel shows a comparison of the lines with those of Orion 


(multiplied by 10, shown in dotted black.). 


in green) is similar to the spatial distribution of different dust 
types identified by Rho et al. (2008; see Figure 3(f)): silicate 
(21 uum dust), carbon, and Al;O; dust regions. These are zones 
of different nucleosynthesis and hence correspondingly very 
different in dust formation and composition. The CO regions 
are in He/O/C and O/C layers (zones 4A and 4B; Sarangi & 
Cherchneff 2013, 2015), and the Ar ejecta are in Si/O layers 
(zone 1B). 


3.2. IFU NIR Spectroscopy 


The positions of the NIRSpec-IFU apertures are shown in 
Figure 3, and representative extracted NIRSpec-IFU spectra are 
shown in Figures 4 and 5. These demonstrate the dominant 
nature of CO emission between 4 and 5 um as well as the 
presence of other emission lines. The strongest ejecta lines are 
from [Ca IV] at 3.2 um, [SiIX] at 3.94 um, [Ca V] at 4.16 wm, 
[Mg IV] at 4.49 um, [Ar VI] at 4.53 um, and [K m] at 4.6 um. 
The fundamental CO features are very bright from 4.2 to 
5.3 um and reveal unique shapes of CO, showing a few tens of 
sinusoidal patterns of rovibrational lines with pseudocontinuum 
underneath. There are many other faint lines that are yet to be 
identified, which is out of the scope of this Letter. There are 
strong lines of [Ar VI] 4.53 jum on the top of the CO feature. 
We also note the wavelength range within the MIRI F770W 
image is dominated by [Ar II] at 6.985 um, which is 100-1000 
times stronger than the continuum (Ennis et al. 2006; Rho et al. 
2008; Smith et al. 2009). 

The strongest detected lines and their properties are 
summarized in Table 1 with derived velocity information. In 
contrast with the previous Spitzer data where multiple velocity 
components were blended together, JWST resolves the knots in 
velocity space. The northern field shows pairs of blueshifted 
emission lines at velocities of —2000 and —1200 km s. The 
southern field shows a single velocity component at 
--1600 km s^'. The observed Doppler velocity structures are 
consistent with the 3D plot (Doppler velocity vs projected 
radius) of ejecta from the Spitzer Infrared Spectrograph (IRS) 
observations (Isensee et al. 2012; DeLaney et al. 2010). 
However, the JWST spectra reveal two reasonably well- 
separated Doppler-shifted groups of ejecta layers at —2000 and 
— 1200 km s^ ' in the northern filaments and a group of knots at 


+1600kms~' in the southern filament. Interestingly, all 
observed elements of Mg, Si, Ar, K, and Ca ejecta show 
similar velocity components. 

A combination of the velocity shift and width indicates 
different nucleosynthetic layers. The two velocity shift 
components at —1200 and —2000 kms! for the NE region 
are likely caused by the different nucleosynthetic layers (see 
discussion by Isensee et al. 2012). However, the velocity width 
may be only a factor of 2—3 times larger than the instrumental 
spectral resolution (~111kms~'), and yet its error is large 
(20960-10096). There may be some separation between these 
layers by a few hundred kilometers per second (corresponding 
to ~5” or ~0.1 pc), and the width may indicate the physical 
size of the layer. For example, [CaIV] layers are separated 
almost by -800kms ! (—0.3pc) between —2000 and 
—1200kms ! layers with each of ~0.1 pe size, while the 
separation of [K m] layers is unclear due to the significant 
errors in the velocity widths. 

We compare the NIRSpec spectra between the north and 
south regions in Figure 5. When we shift the spectra to align 
the features, we see a velocity shift of roughly 3400 kms ' 
between the north and south. The north is blueshifted compared 
with the south and consistent with previous 3D reconstruction 
of the remnant (DeLaney et al. 2010; Milisavljevic & 
Fesen 2013). We also compared the forest of CO features in 
Cas A with those of the Orion nebula. Cas A's CO features are 
bright with prominent individual rovibrational features and 
show a broad underlying component, which we attribute to 
kinematic and instrumental broadening (see Section 4.2 for 
details). JWST spectra of the protoplanetary disk HH 48 also 
show CO features similar to Orion and Cas A but in absorption 
(Sturm et al. 2023). 

We compare the NIRSpec spectra from a few regions. In 
Figure 6 we show spectra extracted from various subregions 
within the IFU fields. The locations of the spectral extraction 
regions are shown in the top panels Figures 6(a)-(d) and the 
color-coded spectra are shown in the panels below. The overall 
intensity and the contrast of the CO lines show variation, and 
variable strength ejecta lines are seen. The northern position 
spectra show more prominent ejecta lines than the southern 
region spectra. 
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Table 1 
Observed Spectral Line Width and Brightness in the Two Representative Spectra 


Region Line Xo Wavelength FWHM Vel. Shift Vel. Width Line Brightness 
(um) (um) (um) (km s7!) (km sl) (10? erg scm? sr!) 
N [Ca IV] 3.2067 3.1862 + 0.0005 0.0059 + 0.0013 —1918 + 46 552 + 122 7.970 + 1.790 
[Ca IV] 3.2067 3.1937 + 0.0005 0.0033 + 0.0010 —1216 + 47 309 + 94 3.948 + 1.196 
N [Si Ix] 3.9357 3.9086 + 0.0019 0.0069 + 0.0047 —2065 + 145 526 + 358 2.669 + 1.839 
[Si Ix] 3.9357 3.9190 + 0.0016 0.0052 + 0.0038 —1151+121 396 + 290 2.123 + 1.563 
N [Ca V] 4.1585 4.1315 + 0.0016 0.0062 + 0.0038 —1947 + 115 447 + 275 2.475 + 1.531 
[Ca V] 4.1585 4.1419 + 0.0013 0.0043 + 0.0032 —1197 + 94 310 + 230 1.709 + 1.261 
N [Mg IV] 4.4866 4.4583 + 0.0019 0.0061 — 0.0049 —1892 + 70 407 + 327 5.230 + 4.109 
[Mg IV] 4.4866 4.4683 + 0.0019 0.0047 + 0.0046 —1223 + 127 314 + 307 3.587 + 3.515 
N [Ar VI] 4.5295 4.4990 + 0.0012 0.0073 + 0.0028 —2020 + 80 483 + 185 11.160 + 4.363 
[Ar VI] 4.5295 4.5105 + 0.0010 0.0049 + 0.0024 —1258 + 1192 324 + 150 7.048 + 3.505 
N [K m] 4.6180 4.5870 + 0.0046 0.0046 + 0.0044 —2013 + 278 298 + 285 3.365 + 4.970 
[K m] 4.6180 4.5985 + 0.0067 0.0070 0.0170 —1266 + 435 454 + 304 2.594 + 6.279 
S UID? s 3.1055 + 0.0001 0.0086 + 0.0003 ses 830 + 10° 2.020 + 0.083 
S [Ca IV] 3.2067 3.2236 + 0.0000 0.0050 + 0.0001 +1581 + 09 468 + 10 2.922 + 0.125 
S UID* 3.7547 + 0.0002 0.0148 + 0.0006 1182 + 48^ 0.861 + 0.052 
S UID* 3.9024 + 0.0001 0.0035 + 0.0002 263 + 15 0.642 + 0.045 
S UID* 4.0570 + 0.0000 0.0050 + 0.0001 365 + 8 4.082 + 0.163 
S CO v= 1-0 E 4.3248 + 0.0002 0.0100 + 0.0004 is 634 + 284° 3.053 + 0.144 
S [Mg IV] 4.4866 4.5109 + 0.0001 0.0063 — 0.0002 1619 + 07 421+ 13 39.080 + 1.955 
S [K m] 4.6180 4.6421 + 0.0003 0.0035 0.0006 1598 + 18 231 + 39 3.021 + 0.707 
Notes. 


* UID indicates an unidentified line. 
> The lines are broader, likely due to a contribution from another faint line(s). 


* The line is broader because of contributions from a few rotational (J) transitions. 


4. Discussion 
4.1. Ejecta and CO Maps 


We have examined the distribution. of Ar, CO, and 
synchrotron emission, shown in Figure 2. Figure 3 shows an 
enlargement toward the bright northern and southern limbs, 
where CO emission can often be seen exterior to the Ar ejecta 
emission. Some examples of regions where enhanced CO 
emission lies at larger radii from the explosion center are 
marked with arrows A, B, and D in Figure 3 and box E in 
Figure 3. We have examined the CO and Ar emission 
morphology and found almost one-to-one correspondence 
between their emission, but CO emission is often shifted to 
appear at larger radii than the Ar emission. We quantified this 
using a radial profile cut across the southern shell (marked as 
Region E in Figure 3). The radial profiles of the cut across the 
southern shell in Figure 7 show that the peak in CO (dotted line 
in green) is located radially outward from that in Ar (dotted line 
in red). Synchrotron emission (in a blue curve) shows a weaker 
correlation with Ar ejecta and CO emission as its morphology 
differs from Ar and CO maps and knots. 

Two scenarios are possible for the observation that the CO 
peak is offset outward from the Ar ejecta. The first scenario is 
that the CO molecules are re-forming behind the reverse shock, 
as suggested by Biscaro & Cherchneff (2014), as the shock 
propagates inwards. The other possibility is simply that Ar 
ejecta has cooled faster than CO, causing a lack of infrared 
emission of the Ar ejecta in the outskirts. In the reverse shock, 
the initial shock velocity is low and the gas density is high, so 


cooling may occur behind the shock wave (Chevalier 1977). 
Since the ejecta are at fairly high temperatures (>10° K), they 
may cool faster than the CO molecules. 

However, the reverse shock moves rapidly as the density of 
expanding gas decreases, so the shock wave soon heats the gas 
to a temperature at which its cooling time is long. Only a few 
percent of the ejected mass may be cooled after the reverse 
shock has passed, as we observe ejecta knots in the far outer, 
high-velocity part of the shell (Fesen 2001). Most ejecta in the 
reverse shock remain at a high temperature, and the lack of 
ejecta emission relative to CO is not due to cooling of ejecta. 

We compare the morphology of 4.5 um (4.3-5.2 jum) CO 
images with the IFU images of [CaIV] at 3.2 um, [Si IX] at 
3.9 um, [MgIV] at 4.5 jum, and [Ar VI] at 4.52 um. Figure 8 
shows separate ejecta distributions of Ca, Si, Mg, and Ar 
emission compared with CO emission in the northern 
NIRSpec-IFU position. The Mg map shows more emission in 
the northern part and different morphology than the others. The 
Si originates interior to the layers than C, O, or Mg. Figure 9 
shows the spatial distribution of the two blueshifted (see 
Table 1) Mg emission components, color coded as blue and red 
in the figure; the blue-colored component stretched north to 
south, and the red-colored component stretches from southeast 
to northwest. Three color maps of [Mg IV], CO, and [Si Ix] 
emission in Figure 9 show contrasting distributions, as the Si 
map differs from the other two. Figure 9 shows the southern 
IFU position, where again the Si map looks different from the 
CO and Mg ejecta emission distributions. The contrasting 
distribution of [Si rX] could be due to the fact that the [Si Ix] 
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Figure 6. (Top) Regions for the extracted spectra shown below are marked on the images of NE IFU FOV for (a) 3.92 um ([Si IX]), (b) 4.70 jum (CO), and of the S 
IFU FOV of (c) 3.92 um ([Si IX]) and (d) 5.04 um (CO). NIRSpec-IFU spectra of Cas A toward a few positions in the northern (middle: (e)) and south (bottom: (f)). 
The positions and marked circles on (a) and (b) and extracted spectra (e) are (R.A., decl.) = (350°8724, 58°8423) in green, (350°8732, 58°8421) in blue, (350°8729, 
58°8419) in red, and (350°8724, 58°8418) in yellow. The positions and marked circles on (c) and (d) and extracted spectra (f) are (R.A., decl.) = (350°8763, 58°7907) 
in magenta ((c), (d)) or orange (f), (350°8764, 58°7904) in cyan, (350°8771, 58°7904) in red, (350°8759, 58°7902) in green, and (350°8760, 58°7900) in blue. 


probably comes from a region of much lower density than the is constant. More extensive IFU coverage of Cas A can help 
lower-ionization lines since the temperature is inversely advance our understanding of the relative locations in ejecta, 
proportional to the density (no) in a shock if the ram pressure CO, and dust. 
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Figure 7. Profiles of surface brightness in Ar ejecta, CO, and synchrotron 
emission along a slit that cuts through the southern shell in an approximately 
radial direction (the slit is marked as box E in Figure 3, bottom panel). The 
position extends from (R.A., decl.) = (3507838, 587789) to (3507836, 587784). 
The CO peaks farther outward than the Ar emission. The vertical lines are 
peaks of Ar (in red) and CO emission (in green), respectively. 


4.2. CO Modeling and Properties 


The Cas A spectra in Figure 4 show two large rounded peaks 
corresponding to the R and P fundamental branches at 4.58 and 
4.83 um, respectively (Rho et al. 2012; see the discussion 
therein). The interval (wavelength difference) between R and P 
branch peaks depends on the temperature. The width between 
them is narrower for a lower temperature (Cami et al. 2010; 
Rho et al. 2012). The spectra show about 70 sawtooth-like CO 
spectral features in Figure 4. The contrast in the CO lines in the 
northern region spectra (red in Figure 5) is higher than that seen 
in the southern spectra (green). The CO spectral features 
largely correspond to specific v= 1—0 vibrational transitions 
from various rotational levels, J. Intrinsic and instrumental 
broadening of the lines blends then with weaker lines from 
v — 2-1 and higher vibrational states and with each other for 
higher J values in the R branch. 

The CO fundamental lines have a similarity to the JWST 
spectra seen in the Orion nebula (blue in Figure 5), indicating 
similar individual CO features between Cas A and Orion. 
However, Cas A has much larger kinematic broadening, which 
changes the spectral appearance. The wings of the lines, and 
weaker lines, blend together to creating a pseudocontinuum. 
Recent JWST observations of HH 211 detect CO features (Ray 
et al. 2023), somewhat similar to those of Cas A, but its 
pseudocontinuum is not apparent like in Cas A. The projected 
dynamic velocity in HH 211 is 80-100 km s '!, which can be 
the upper limit of the CO line width, and the line width is likely 
much less than that in Cas A. 

We performed LTE model fitting based on Cami et al. (2010) 
and the CO molecular line data from Goorvitch (1994) to 
determine the physical conditions of CO molecules. This 
approach is similar to that used by Rho et al. (2012), Banerjee 
et al. (2016), and Rho et al. (20182, 2021). We have quantified 
the goodness of fit (x) for the models using the IDL routine 
MPFIT and the Python routines of MPFIT and LMFIT for the 
best-fitting parameters (Markwardt 2009). Before we compared 
with the CO LTE models, we shifted the velocity of the CO 
features to match the red- and blueshifted fine-structure line of 
[Ar VI] at 4.53 um (we used the brighter component among the 
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two velocity components). The northern and southern knots are 
blue- and redshifted, respectively. This shift allows the CO 
models to match the individual rovibrational lines. 

Figure 10 shows the best fit of an LTE model with a 
temperature of 1078 (900-1250) K, a velocity width (FWHM) 
of 237 (160—350) kms’, and a column density of 3.8 
(1-20) x 10" cm ?. The statistical errors are ~2%, but the 
errors include systematic errors by using different routines and 
methods. The intensities of the CO lines are well reproduced by 
the model, but the contrast in the lines is lower at higher J in 
the R branch, and the level of the pseudocontinuum falls well 
below the data at high J (the extreme long and short 
wavelengths of the emission). The spectral resolving power 
of the IFU is ~1000 at these wavelengths, which translates to a 
velocity resolution of ~111kms~'. However wider Gaussian 
profiles with a velocity width of 237 (160—350) km s! were 
able to reproduce the observed contrast of the CO lines and the 
blended emission underneath, as shown in Figure 10. Isolated 
far-IR CO rotational lines have shown similar velocity widths 
(Wallstróm et al. 2013). After we used line-width broadening, 
no extra smoothing was needed. The fundamental band is a 
combination of hundreds of individual lines, while the far-IR 
CO rotational lines are well isolated at each wavelength. The 
best-fit model (red in Figure 10) with a velocity width of 
237 km s! has too little contrast in the lines compared to the 
observed spectrum (green in Figure 10) in the R-branch CO 
features («4.67 jum), which suggests that a lower velocity 
width may be appropriate. 

The ejecta and molecular knots have different density 
regions within them. The velocity width of ~237 km s^! 
indicates that the CO gas is the post-shocked region of the 
reverse shock. The dense knots are compressed and fragmented 
by instabilities arising from the shock passage (Mellema et al. 
2002; Raga et al. 2007; Silvia et al. 2010). The turbulence in 
dense knots is small, less than the turbulent velocity of 
10 kms ' in the warm gas (Korpi et al. 1999) 

The observed velocity width of CO is the internal shock 
velocity that results from varying CO density and geometry 
(including the scales unresolved) that lead to variation in the 
shock speed and the direction with respect to the line of sight. 
The shock velocity in CO gas (Vco) is inversely proportional to 
the square root of the density ratio between CO (nco) and 
ejecta (ejecta) and proportional to the reverse shock velocity 
(V) (Vco = Vis X 4/Mejecta/Mco) owing to the conservation of 
energy (Biscaro & Cherchneff 2014; Klein et al. 2003). The 
reverse shock velocity V,, relative to the ejecta has beenesti- 
mated to be 2000-3000 km s^! using optical data (Morse et al. 
2004; Biscaro & Cherchneff 2014). 

The reverse shock velocity in the ejecta can also be estimated 
from an X-ray gas temperature of 2.6 keV (~3x10’ K; Hwang 
et al. 2004). The X-ray gas temperature T* (= 3/16 u/k V2 ~ 
11 Và: references in Rho et al. 2021) infers V,, of 1700 km s! 
(Equation 12—24 in Spitzer 1978) where u is the mean mass per 
particle and k is the Boltzmann constant). The density of the 
ejecta is 10^? cm? (DeLaney et al. 2010), and the den- sity of 
CO gas is approximately 10977 cm ? (Biscaro & Cherchneff 
2014; Docenko & Sunyaev 2010). The critical density of CO 
(Nerit = Au/^, Where the Einstein A is the radiative transition, 
and y is the collision coefficient between upper and lower 
levels) is high (> 10° cm? ) for a temperature of 1000 K and 
an example of the upper and lower levels of v = 21 — 20'. 
The observed velocity width of 237 (160—350) km s! is 
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Figure 8. Cas A band images of the northern IFU field: (a) [Ca IV] at 3.2 jum, (b) [Si IX] at 3.9 jum, (c) [Mg IV] at 4.5 jum, (d) [Ar VI] at 4.52 jum, (e) CO 4.79 
(4.788—4.798) um, and (f) synchrotron-subtracted CO broad band (from F444W) images. The latter captures the brightest CO band. The CO emission appears similar 
to other ejecta maps. 
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Figure 9. Cas A band images of the northern IFU field. (Left) Red and blue colors denote the two different kinematic components of [Ca IV] 3.2 um in this field. 
(Middle) Three-color images showing [Mg IV] (in blue), CO (green), and [Si IX] (red). (Right) Three-color images of the southern IFU field: [Mg IV] (in blue), CO 
(green), and [Si IX] (red). CO and [Mg IV] show similar morphology, while [Si IX] morphology is different. 
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Figure 10. An LTE model (in red) using a temperature (T) of 1078 K and a velocity width of 200 km s^! is superposed on the JWST southern region spectrum of Cas 
A (in green). The sawtooth-like patterns are reproduced, but the relative strengths differ slightly. In the expanded view in the expanded view of the spectra (Right), the 


peaks are from CO v 


approximately the shock velocity inside the CO clumps of 
170—300 km s~' (= 1700-3000 x J6 x 105/(5 x 109). 

The LTE model to reproduce the JWST spectra of Cas A 
only uses 12C!60 without other CO isotopes. We simulated the 
intensity profile for optically thin and thick cases and found 
that the CO emission in Cas A was close to the optically thin 
case as it was found in the low-resolution AKARI spectra (Rho 
et al. 2012). The CO column density of 3.8 (1-20) x 10" cm ? 
is comparable to those obtained from the CO excitation 
diagram (Wallstróm et al. 2013). The high column density is 
consistent with the idea that CO molecules are re-formed 
behind the reverse shock (Wallstróm et al. 2013; Biscaro & 
Cherchneff 2014). The corresponding CO mass of the knot in 
the southern field is ~2.3 x 10 Mo per pixel (0/1) at a 
distance of 3.4 kpc. Eventually, summing up all CO emissions 
from F444W emission and estimating the total CO mass from 
Cas A will be an interesting value. For that, JWST’s advanced 
pipeline processing of individual frames is required, including 
removing hot pixels due to cosmic rays and various artifacts 
related to star subtraction, which is out of the scope of this 
Letter. The dip between the R and P branches (around 4.65 jum) 
from the LTE model is slightly off from the JWST spectra that 
were shifted based on ejecta lines. When we made a further 
shift of —0.014 um (a total of 0.03 um), the dip in the model 
was consistent with the JWST spectra of Cas A. 

We notice a difference between the CO model and the JWST 
spectra at shorter and longer wavelengths as shown in 
Figures 10 and 11, which is much larger than the synchrotron 
emission (—5 MJysr !) The synchrotron emission is an 
interpolation from radio fluxes, and detailed descriptions can 
be found in Rho et al. (2003, 2008, 2012). The emission could 
be due to dust emission. Therefore, we added a continuum 
component using carbon dust to the CO component and fit the 
spectrum. The results show the continuum at 4.8—5.2 um has 
improved (the reduced x improved a factor 10), but the 
continuum fit did not improve the residuals between 4.25 
and 4.4 um. 

An interesting feature of the Cas A CO spectra is the band 
heads that appear at ~4.3 um. These band heads are a rare 
detection. Starting from J — 0 at ~4.6 um, the wavelengths of 
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1-0 P(4) through P(32) transitions. Blends with weaker v = 2-1 and higher vibrational states modulate the apparent contrast of the lines. 
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Figure 11. A two-temperature LTE CO model (in red) is superposed on the 
JWST southern spectrum of Cas A (in green). Here we have added a hotter CO 
component with T= 4800 K (yellow) and AV —210 kms ' to the first- 
temperature component (blue). The total fit is shown in red. The band heads 
(v = 1-0, v = 2-1, and v = 3-2) are marked. The contribution of synchrotron 
emission is negligible about ~50 MJy sr !, and wavelength dependence is 


small (the slope is 0.67). 


successive R branch transitions get progressively shorter as J 
increases. However, the lines become increasingly closely 
spaced, and after J 90 the lines begin increasing in 
wavelength. The unresolvable pileup in wavelength of lines 
near J ~ 90 creates the band heads. Cas A spectra clearly show 
these band heads for at least three different vibrational 
transitions (v — 1-0 at 4.3 um, v —2-1 at 4.35 um, and 
y — 3-2 at 4.45 jum), marked in Figure 6(e). The laboratory 
band-head wavelengths are v= 1-0 at 4.295 um, v — 2-1 at 
4.350 um, v = 3-2 at 4.295 um, and v= 4-3 at 4.465 um 
(Geballe et al. 2007). The observed band heads are shifted by 
about —2000 km s^! and +1000 km s^! for the northern and 
southern region spectra, respectively, as we have previously 
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noted. The excitation energy of the v = 0 J = 90 level is about 
~22,000 K, which implies that the temperature must be high 
for even a small fraction of the CO to be excited to such high v 
and J states. Therefore, the presence of such band heads may 
indicate multitemperature CO gas, or it may require a non-LTE 
model. We attempted to fit a two-temperature CO LTE model 
(with a second temperature component of 4800 K), shown in 
Figure 11 although the fit did not converge efficiently due to 
many small sinusoidal patterns and diverse combinations 
of two CO models. However, it demonstrates that the high- 
temperature CO model can reproduce the CO band heads at 
4.3, 4.35, and 4.45 m, indicating multiple temperature 
components of CO gas likely exist. 

We also analyzed the CO excitation with the non-LTE 
radiative transfer code RADEX (van der Tak et al. 2007) as it 
was done by Wallstróm et al. (2013). The more sophisticated 
methods retain the assumption of a local excitation but solve 
for the balance of excitation and de-excitation rates from and to 
a given state, the so-called statistical equilibrium. The RADEX 
program was adapted for mid-infrared (MIR) and NIR CO 
emission (Wallstróm et al. 2013). While the model reproduced 
low-resolution CO AKARI spectra using temperatures of 400 
and 2000 K, a velocity width of 300 km s |, and Nco of 
4 x 10!” cm”, it could not reproduce the CO sawtooth-like 
patterns seen in the JWST spectra. The contrast in the CO lines 
was much larger than in the JWST spectra (more than a factor 
of 5 stronger) when using a velocity of ~200-600 km s, and 
increasing the velocity did not help reproduce the JWST 
spectra. The RADEX model produces asymmetrically oscillating 
CO patterns, which differ from the observed JWST spectra. 
The model, in particular, lacks CO emission below 4.5 jum and 
above 4.9 jum, which could be due to lack of H-atom collision 
partners or other missing collision rates. Future work on non- 
LTE modeling is encouraged to reproduce the JWST spectra of 
Cas A, which may need to include a large Doppler velocity 
broadening and various collision rates, including H-atom 
collision partners. 

Non-LTE CO modeling can be preformed using the 
hydrodynamical radiation code (HYDRA), a hydrodynamic 
SN evolution code, including the nuclear networks and ejecta 
mixing, and Radiation code, accounting for the atomic-level 
population and the equations of state (see Hoflich et al. 1993; 
Rho et al. 2021). The y-rays and free electrons at a certain 
temperature can ionize CO, which may affect the timescale of 
the CO and CO* formation and the ratio of CO* and CO. The 
time-dependent rate equations are C +O — CO and Ct + 
O — CO* (see also Liljegren et al. 2023; Ono et al. 2024). 
Nonthermal excitation can enhance CO*, which adds a 
blueshifted emission component. NIR and MIR spectral 
signatures from CO only and from CO and CO* with and 
without nonthermal excitation for Type IC SNe are described in 
the appendix of Rho et al. (2021). Such complex HYDRA 
modeling for the progenitor Type IIB SN is encouraged and 
can be applied to the JWST spectra of Cas A. 


5. Conclusion 


In this Letter we have presented and analyzed JWST 
observations of the SNR Cas A to investigate the physical 
conditions for molecular CO formation and destruction in SN 
ejecta. Our conclusions can be summarized as follows: 
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. NIRCam (F356W and F444W) and MIRI (F770W) 
images map synchrotron radiation, Ar-rich ejecta, and 
CO in great detail and reveal complex structures on the 
scale of individual ejecta knots. The CO map strongly 
correlates with knotty ejecta structures. CO is brightest in 
the northern and southern portions of the main shell and 
weakest in the east where ejecta emission dominates. 

2. Complex fundamental bands of CO features dominate the 
NIRSpec-IFU spectra, with ~70 individually resolvable 
rovibrational lines. These lines are similar to those seen in 
the Orion nebula but with instrumental and velocity 
broadening in Cas A that blends lines to produce a 
pseudocontinuum. 

3. LTE modeling of the CO spectra (R and P branches) 
suggests an optically thin CO gas with a temperature of 
~1080 K, N(CO) = 4 x 10" cm ?, v 2200 kms |, and 
Mco —2 x 10^" Mo per 0/1 pixel. 

4. We discover band heads for the v — 1—0, 2-1, and 3-2 

rovibrational transitions at 4.3, 4.35, and 4.4 um. The high 

rotational levels (J ~ 90) that need to be populated to form 
the band-head emission requires a high-temperature 
component in modeling the emission. When we fit the 
spectrum including a second CO component, the second 
temperature is ~4800K, and the fitting favors a lower 

first-temperature component (ie. lower than 1080 K), 

indicating the presence of multiple temperatures. 
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